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ABSTRACT 

Faster-than-light or superluminal motion was originally predicted as a relativistic illu- 
sion of ballistic moving ejecta, and confirmed in a few tens of sources observationally. 
However, the recent results of the long-term multi-epoch observations of quasars, ac- 
tive galaxies, tracing the structure further along the jets and following the motion of 
individual features for longer time, rise questions that are difficult to understand by 
the standard ballistic model. I.e., the ejecta are aligned with the local jet direction, 
instead of the core; and within individual jets apparently inward-moving features are 
observed. Here we show that these unexpected phenomena, although only a small frac- 
tion among large samples, indicate the existence of non-ballistic jet motion, in which 
a continuous jet produces a discrete hot spot. And the precession of such a hot spot in 
the plane of the sky appears superluminal. Therefore, an unified and simple interpreta- 
tion to the new results is obtained, which can be further tested through its predictions 
on the evolution of ejecta. The study is of importance in the understanding of the 
nature of superluminal motion, the interaction of jets and surrounding materials, as 
well as the common physics underlying quasars and microquasars. 

Key words: (galaxies:) quasars:general-(stars):binaries:general 



1 INTRODUCTION 

The appearance on the sky of relativistically moving out 
flow, which expands at speed greater than the speed of 
light, has been predicted by Rees (1966), five years in ad- 
vance of the discove r y of superluminal motion (|Rees 19661 ; 
IWhitnev et al. 197ll ; ICohen et al. 197ll ). In units of the 
speed of light, c, the apparent speed j3 a = v/c of separa- 
tion is the ratio of the difference in observed positions to 
the observed time-interval. With the blob travelling at an- 
gle 9 with respect to the line of sight (LOS), j3 a is given 
by 



Pa 



P P sin ( 



1 — f3 p cos 8 



(1) 



When v ~ c, and 6 at certain range, the apparent speed, 
f5 a , can be larger than the speed of light. In fact Eq. JTJ and 
Fig. la correspond to a simple ballistic motion of ejecta, with 
a fixed jet axis, which has been applied to the superluminal 
motion of quasar like 3C279. 

To interpret the S-shape patte rn of SS433, a modifie d 
ballistic model has been proposed (|Abell fc Margon 19791 ), 
in which the direction of the jet axis is variable, and denoted 
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by a, b, and c in Fig. lb. In this model, features b and b' 
in Fig. lb, are the same ejection evolving at different times, 
which should always move along the jet axis o — b. As the jet 
axis precesses to another direction, i.e., o — c, another fea- 
ture ejects out. Similarly, the trajectory of it should always 
be fixed at that direction. Assuming that only one feature is 
ejected at each direction, then the connection of these fea- 
tures ejected at different directions like a, b and c forms a 
one sided spiral pattern, as in Fig. lb. And considering the 
receding features at opposite direction to b and c a S-shape 
pattern, like SS433, is explained. 

Generally, the ballistic model, with jet axis fixed or not, 
implies that each individual feature must move along the jet 
axis away from the core. Thus, the observed inward motion 
of features can be interpreted by newly emerging jet features 
ejecting from the core at the right position and time, which 
mimics the decrease of the apparent separation of the core 
and other features. 

Therefore, the ballistic model implies that the whole 
curved pattern would move further away from the core at 
the next step of evolution due to that each individual feature 
evolves linearly as in Eq. {TJ. However, this is not consistent 
with the observation that the jet flow apparently occur s 
along preexisting bent channels (lKellermann et al. 20041 ). 
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Figure 1. A schematic illustrating of three models of super liminal 
motion. The first one is the pure ballistic model with a fixed jet 
axis, shown at the left side. The second one is the model of ballistic 
plus precession, represented by the solid spiral connecting the 
features, a, b, and c. And the last one is the non-ballistic model, 
denoted by the dashed ellipse through feature b, which precesses 
with a constant distance, Rd, to the core. 



instead of away from the core as a whole. Therefore, it is 
conceivable to consider other possibilities of jet motion. 



2 THE NON-BALLISTIC MODEL 

Theoretically, a hot spot can be generated by the interaction 
of a continuous jet with the surrounding material. The pro- 
cess of continuous ejection producing "discrete hot spots" 
can be simply extracted from Dermer (1999). The energy 
of the injected nonthermal particles comes at the expense 
of the directed bulk kinetic energy of the fluid. Assuming 
that the system produces an outflow with total energy Eo 
and initial bulk Lorentz factor Fq. Since most of the energy 
of the flow is bound up in the kinetic energy of baryons, 
assumed to be protons, t hen Eq = ToN t m v (? , with Nt the 
total number of protons jDermer 19991 ). 

Through the dynamics of the blast wave, which deceler- 
ates by sweeping up material from the surrounding medium, 
the deceleration radius, which can be consid ered as the dis - 
tance of a hot spot to the core, is obtained as (jDermer 19991 ), 



R 



(3 - C)Eo 



,1/3 



4itf b n r 2 m p c 21 ' (2) 

where the density of external medium can be parameterized 
by the expression n ex t{R) = na(R/ Rd)~ C ' ■ In the simplest 
case, one finds that T(R) Si T for R < R d , and T(R) oc R~ 3 
for R 3> Rd, where g — 3 — ( and g = (3 — C)/2 in the 
adiabatic and radiative regimes, respectively. The term 
represents the fraction of the full sky into which the explo- 
sion energy is ejected. The decelera tion radius, Rd, occurs 
at the observing time (jDermer 19991 ), 



t d = r (l - A) cos0)(l + z)R d /(cT ), 



(3) 



where (3q = yl — r ( 7 2 and z is the redshift and 6 is the 
misalignment angle between the jet and LOS. At a given 
frequency and 6, there is a peak in flux density at the time 



td of Eq. <(3j , in other words, the brightest emission region 
during its whole evolution time appears at td and distance 
Rd to the core. 

The light curves for the model synchrotron flare at X, 
gamma, and radio frequencies, both along the jet axis (thick 
curves), and at 20deg to the jet axis (thin curves) are shown 
in Fig. 2 of Dermer (1999). Obviously, each peak in the light- 
curve corresponds to a hot spot during its whole evolution 
time. 

The radial size of such a hot spot can be defined as: the 
length corresponding to the decline of the flux density from 
the peak for one order of magnitude (for X-ray emission), or 
for a factor of five below the peak (for radio emission) . Thus, 
the radial length of a hot spot compared with its distance to 
the core is SRd/Rd = St/t, which is between 10" 1 to 10~ 3 
in Fig. 2 of Dermer (1999). 

When the jet axis precesses to another direction, the 
process of ejection and deceleration repeats. Hence another 
hot spot is produced, and its distance to the core is still R 
in the case of an isotropic distribution of the material. Thus 
different hot spots are generated by the same continuous jet, 
interacting with matter at different directions, which mimics 
one hot spot moving in space continuously. 

Contrarily, for the ballistic model hot spots or fea- 
tures are produced by discrete ejection events. The dis- 
tance of each feature to the core is determined by Eq. (JTJ, 
which varies with time. Thus simply replacing the ra- 
dial v elocity of a feature, v, of the ballistic model o f 
SS433 (|Abell fc Margon 19791 ; [Hiellming fc Johnston 198lh . 
by the distance of it to the core, Rd, the precession of a hot 
spot under the non-ballistic model is obtained as, 

Rx = Rd [sin A sin i cos r\ + cos A cos i] , 
R y = Rd [sin A sin rj\ , 

Rz = Rd [cos A sin i ~ sin A cos i cos rj\ , (4) 

where R x , R y , and R z are the components of Rd in the coor- 
dinate system x — y — z. The :r-axis is towards the observer, 
rotating around the a;-axis for angle £, so that the new y-axis 
(AS) will point north, and the new z-axis (Act) will point 
east. As shown in Fig. 1, A is the opening angle of the preces- 
sion cone, i is the inclination angle between the jet rotation 
axis and LOS, and the precession phase is r\ = tit + 770 (770 is 
the initial phase). The displacement of a feature to the core 
can be described with respect to AS and Act as, 



R a = R y sin £ + R z cos £ , 
Rs = R y cos £ — R z sin £ . 



(5) 



Differentiating Eq. the velocity projected to the plane 
of the sky is given by, 



v a = R y sin £ + Rz cos £ . 
vs = R y cos £ — Rz sin £ . 



(6) 



The time taken by a spot to precess for tangent distance 
of the size of a hot spot, SRd, is St = Srj/Q, where Sr) — 
SRd/(R d sin A). The cooling time of a hot spot, St co , corre- 
sponding to, i.e., the time taken for the radio peak to decline 
for a factor of five, can be inferred from Fig. 2 of Dermer 
(1999). If St < St co , then the hot spot appears as a filament, 
otherwise, it appears as a hot spot. Obviously the precession 
can make a filament bent or twisted. 
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3 THE APPARENT LUMINOSITY 

With the intrinsic luminosity is Lo , the apparent luminosity 
of a jet feature is given by, 



L S" 



(7) 



where 8 — 7 i r 1 (l — (5b cos^) -1 , and 75 is the Lorentz factor, 
given by 76 = (1 — /3 6 2 )~ 1/ ' 2 . And n depends on the geometry 
and spectral index, which is typically in the range between 
2 and 3. The angle 9, denoting the misalignment angle be- 
tween the LOS and the jet axis, changes as the precession 
of the jet axis, 



= cos A cos i + sin A sin i cos 77 . 



(8) 



As 8 varies, a feature appears moving away from the core 
or close to the core at different precession phase, as shown 
in Fig. 1. Hence, the apparent luminosity varies with the 
precession. Obviously, the maximum luminosity occurs when 

e ~ 0. 

Whether a feature is detectable or not is dependent of 
the misalignment angle between the jet axis (correspond- 
ing to the hot spot) and the LOS, 8, as given by Eq. (|8]). 
The larger the angle, 8, the weaker the Doppler boost- 
ing effect. The jet emission becomes undetectable when 
\8\ exceeds certain limit. This explains the disappearance 
of features at large angular disp lacement, i.e., for GRS 
1915+105 l|Miller- Jones et al. 20071 ). 

By the new model the apparent velocity of a hot spot, 
/3 a , can be represented as, 



Pa oc 



SIR sin A 



(9) 



which is dependent of parameters such as: R, SI, A, i, and 77. 
Thus, the corresponding pattern speed, j p , if inferred from 
Eq. {TJ, is not directly related with the bulk speed, 7b . This 
explains the observation of both the sl ow quasars with 7„ <C 
7i, and the fast quasars with 7 P ~ 76 l)Cohen et al. 2007 



4 INWARD MOTION AND BENT 
TRAJECTORY 

The apparent inwar d motion for 

quasars dKellermann et al. 20041 ) can be explained in 
the context of non-ballistic model. For some samples, the 
LOS is as the left side of Fig. 2, denoted by LOS\, in which 
ejecta are moving out of the core (A2 and B2 are projected 
to 1 and 2, respectively); for other ones, the observer's view 
is as LOS2, at the right side of Fig. 2, where the ejecta are 
moving towards the core (B2 and C2 are projected to 3 and 
4, respectively). Under such a circumstance, spot C2 would 
precess to D, which appears moving toward to the core, by 
projecting to D' on the plane of sky. 

Beside a qualitative interpretation of apparently inward 
motion, as in Fig. 3, the observation of quasar 0119+041 can 
also be fitted by the new model. Assuming the parameters 
SI = 43.0deg/yr, R d = 2.9mas, i = ll.Odeg, A = 3.7deg, 
rj — 225.0deg, and £ = 4.3deg, and putting them into 
Eq. ([S}, Fig. 3 can be obtained. Obviously, the curve given 
by the new model fits the observations (dots) better than 
the dashed line of the ballistic model, shown in Fig. 3. The 




Figure 2. A schematic illustrating of inward motion under the 
non-ballistic model. The precession of jet axis results in motion of 
the brightest feature in the elliptic trajectory. For some sources, 
the observer's view angle is as the left side which appears that 
ejecta are moving out of the core (A2 and B2 are projected to 1 
and 2 respectively); whereas for other sources, it is as the right 
side in which ejecta seem to move inwardly to core (B2 and C2 
are projected to 3 and 4 respectively). The next evolution would 
be that C2 moves to D which is projected to D' which is even 
closer to the core. 




Figure 3. Fitting of 0119+041. The observational inward-moving 
features of 0119+041 are represented by dots, which are fitted by 
the curve obtained by putting the parameters of 0119+041 in the 
text into Eq. J5]l. The dashed line is the fitting by the ballistic 
model. 



precession velocity of the jet axis, SI, can originate from ei- 
ther the Lense-Thirring effect, or the torque exer ted by the 

companion object in a binary syste m (k atz 1997^ 

The observed bent trajectories (|Kellermann et al. 2004T) 
can also be explained by the new model. According to 
Eq. ©, the whole trajectory of a hot spot on the plane 
of the sky in one precession period is an ellipse. And part of 
this elliptic trajectory appears as a bent trajectory, which 
should be nonlinear and non-radial. This explains why the 
ejecta are aligned with the local jet direction (as 1-2, or 3-4), 
instead of the core (o-2, or 3-o), as shown in Fig. 2. 
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Figure 5. The multi-ejection of GRS 1915+105 can be fitted by 
the parameters of GRS 1915+105 in Table 1, in which the three 
sub-hot spots satisfy A fs i, so that they appear as ejected from 
the core. 



Figure 4. A schematic illustrating of multi-sub hot spot scenario. 
Each one has its own misalignment angle with respect to the jet 
axis, i.e., A2 for sub- hot spot 2, and its distance to the core, R^, 
can also be different. In this figure, sub-hot spot 1 and 3 appear to 
move away from the core, denoted by o, whereas, spot 2 appears 
to move inwardly to the core. 



5 MULTI-EJECTION 

The creation of multi-components in pulse profiles of a pu l- 
sar may result in a nested co ne structure (|Gil et al. 19931) . 
or a patchy beam structure l|Lvne fc Manchester 19881 ). If 
the continuous jet producing hot spots in quasars and mi- 
croquasars has a similar clumps beam, then the interaction 
of such a beam (or structured jet) with the interstellar ma- 
terial can reproduce multi-components either. In such case, 
each component can have its own misalignment angle with 
the axis of the precession cone, and its distance to the core 
can also be different. 

It is conceivable that each component has approxi- 
mately the same precession velocity around the axis of pre- 
cession cone. So that the precession period of the major 
beam around the axis of precession cone and the pattern of 
the major beam are unchanged during jet precession. The 
precession of such a jet mimics a multi-component (or sub- 
spot) moving in space as in Fig. 4. 

By assuming the parameters of Table 1, the multi- 
ejection of microquasar GRS 1915+105, can be fitted as 
in Fig. 5. As the beam precesses through LOS, only 
these sub-hot spots, with the opening angle of preces- 
sion cone, Xt, equal the inclination angle, i, correspond 
to a zero angular displacement to the core, which appear 
to be ejected from the core. Those with 7^ i, pre- 
cess through LOS with nonzero angular displacement and 



Table 1. The fitting parameters of and GRS 1915+105. 



Q 


i 


c 


Ai 


A 2 


A3 


0.283 


52.9 


285.0 


52.9 


52.9 


52.9 


ill 


>I2 


V3 


Ri 


R2 


Rs 


0.0 


-1.2 


-2.6 


5350.0 


5350.0 


6500.0 



Q is in unit deg/day. R%, R2, and R3 are the distances of the three 
sub-hot spots to the core in unit of mas. All other parameters 
are in unit of degree. The receding feature of GRS 1915+105 is 
obtained by assuming R r = 2675. Omas, and r) r = 5.5deg. 

hence don't appear to be ejected from the core, which ex- 
plains the phenomena i.e., on 07 36+017, 0735+178 and 
1219+285 (|Kellermann et al. 2004 ). 



6 PRECESSION OF A BENT JET 

The source 0735+178 is a bizarre object, since all strange 
phenomena occur on it: inward motion (square dots), not 
from the core (all dots) and ejecta aligning with the local jet 
direction, instead of the core (all dots), as shown in Fig. 6. 
Moreover, it contains both outward motion (dots in filled 
and hollow circles) and inward motion (square dots). 

Under the new model, the jet motion of 0735+178 re- 
sults in the precession of a bent jet, which has three hot 
spots instead of one in 0119+041. The three hot spots can 
be seen as Ci, B2, and A3 in Fig. 2, locating at different 
radial distances, Ri, R2, and R3 respectively. The hot spots 
corresponding to Ri and R2 are not bent much, when the 
bulk velocities (radial), Vi/c and V2/C are much larger than 
the precession induced tangent velocities, /3i and /?2, given 
by Eq. ©. As the distance of the plasmoid to the core be- 
comes greater than R2, its velocity drops considerably, so 
that vzjc is less (or much less) than that of precession in- 
duced velocity, say, (3$. Thus the hot spot A3, appears bent 
significantly. The image of 0548+ 165 may provide an ex am- 
ple of such significantly bent jet (jMantovani et al. 19981 ). 

Under the condition of an isotropic distribution of ma- 
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Table 2. The fitting parameters of 0735+178 



n 


% 


i 


Ai 


A 2 


A 3 


19.6 


12.6 


50.9 


9.2 


15.5 


6.9 


m 


m 


m 


Ri 


R2 


R 3 


-154.7 


-166.2 


-17.2 


4.0 


8.4 


13.1 



Q is in unit of deg/yr. Ri is in unit of mas, and others are in unit 
of degree. 

terial, and a continuous precessing jet (period of precession 
is constant), a bent jet would precess in space with an un- 
changed shape, as shown in the three curves through A2, B2 
and C2 respectively in Fig. 2. 

Consequently, the strange morphology of 0735+178 can 
be explained. The three hot spots along the bent jet can be 
treated as three individual straight jets, with different mis- 
alignment angles to the axis of precession cone, and initial 
precession phases, as the dashed lines in Fig. 2. In such case 
the jet motion of 0735+178 can be fitted as Fig. 6 through 
the parameters in Table 2. 

Obviously the motion of ejecta of 0735+178 is more 
complicated than that of 0119+041, because the former con- 
tains three hot spots along a bent jet, whereas the latter 
contains one. 

From the context of the non-ballistic model, there 
are two main discrepancies between GRS 1915+105 and 
0735+178, firstly the former satisfies, A m i, whereas the 
latter satisfies A / i; secondly 6 can only be observed in 
a small range (a few degrees) in the former, whereas it is 
observed in much larger range in the latter. 

The first discrepancy explains why the features of GRS 
1915+105 appear ejecting out from the core, whereas that of 
0735+178 don't. And the second discrepancy explains why 
the trajectories of GRS 1915+105 appear linear and that 
of 0735+178 appear non-linear. In fact a curved trajectory 
in GRS 1915+105 is inevitable, if the ejecta of it could be 
observed at large 6. 

Notice that the fitting parameters of 0119+041 and that 
of Table 1 and Table 2 for GRS 1915+105 and 0735+178 
respectively are not unique. Other numerical solutions to 
them cannot be excluded. 



7 OUTWARD VS INWARD 

The new model can account for the very large numbers of 
observed outward motions, as opposed to the inward ones. 

As shown in Fig. 7, the dash profile corresponds to the 
apparent luminosity, L, obtained by putting of Eq. ([8]| into 
Eq. (0. Since 9 is a function of time (due to 77), L is depen- 
dent of t either. In fact the profile of L(t), or the light-curve 
(L-C), can only be obtained in the case that an emission 
beam with negligible opening angle precesses through LOS, 
in other words, the convolution of L(t) with a S(t) function. 

Generally, the L-C appears as the convolution of L(t) 
with the emission beam, P(t), which has certain opening 
angle. In such case, the L-C is given by the convolution, 

F(t) = J" L(r)P(t-r)dr, (10) 

where P(t) can be both homogeneous or nonhomogeneous. 
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Figure 6. Fitting of 0735+178. The three hot spots along the 
bent jet can be treated as three individual straight jets, precess- 
ing with the same velocity around the axis of precession cone, 
but with different misalignment angles with respect to the axis of 
precession cone, and different initial precession phases. The tra- 
jectories of ejecta are fitted by three curves via the parameters of 
0735+178 in Table 2. 



Here its shape is assumed to be a rapid rise and slow de- 
crease, which can be denoted by i.e., the triangle at the 
bottom of Fig. 7. By t a = — oo and tb = ti, the L-C of an 
inward feature is obtained. Similarly, by setting t a = ti and 
tb = oo, the L-C of the outward feature can be obtained 
either. 

For simplicity, the core emission is assumed to be a con- 
stant, as the horizon line with a hight of unit shown in Fig. 7. 
However the jet emission can be attributed to the emission 
of the core, when the jet precesses to the small 8(t) region 
within the core area. The core area is between the time inter- 
val, ti to <2, as shown in the slash rectangular. Thus the ap- 
parent luminosity of the core area becomes, F(t) + 1, where 
F(t) is given by setting t a = t\ and tb = ti in Eq. (|10[l . 

By Eq. (|10p . the apparent luminosity of the core, the 
inward and the outward features are obtained. The appar- 
ent luminosity of an inward feature increases gradually. It 
disappears in the core area (region 1), at its maximum F(t), 
as shown in the left hand side curve of Fig. 7. The region 1 
is the common area of the core and the inward jet, in which 
the decrease of the jet apparent luminosity occurs during the 
increase of that of the core, which makes this region difficult 
to distinguish from the core. 

Contrarily, an outward feature starts in the region 2, 
which is beyond the true core area. The jet emission rises 
rapidly, after reaching its maximum F(t), it decreases grad- 
ually. 

Therefore, compare with the inward features, the out- 
ward ones have an additional region 2, which is favorable 
to observe. This may explain the phenomenon of very large 
numbers of observed outward motions, as opposed to the 
inward ones. 

The radio L-C of the SE component of GRS 1915+105 
at 8.46GHz (Miller- Jones et al. 2007), as shown in the small 
box (A, B and C) is consistent with the outward curve (A' , 
B', and C) in Fig. 7. 
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Figure 7. The dash profile corresponds to the apparent luminos- 
ity, L(t). The convolution of L(t) with the emission beam, -P(t), 
represented by the triangle at the bottom gives following results: 
the apparent luminosity the core, the inward and the outward 
features, as shown by the top dashed curve, the bottom left hand 
side curve, and the bottom right hand side curve respectively. 
The hight of the light gray rectangular at the bottom represents 
the limit of telescope sensitivity. The core area is denoted by the 
slash rectangular between t\ and t2- The small box at the top 
left is the radio L-C of the SE component of GRS 1915+105 at 
8.46GHz. 



8 DISCUSSIONS AND PREDICTIONS 

The non-ballistic model explains superluminal motion by the 
apparent transverse motion of hot spots, resulted in a con- 
tinuous jet interacting with material at different locations 
and directions at approximately the same distance to the 
core. Obviously, the scenario doesn't mean the transporta- 
tion of energy and moment from one point to another with 
the apparent speed. 

The appearance on the sky of precessing jets is affected 
by time delays, since the image observed is composed of 
photons coming from various parts of the jet that happen 
to have the same arrival times. However, in the fitting of jet 
motion of quasars and microquasars using the non-ballistic 
model of this paper, the time delay is neglected. The reason 
is as follows. 

The time delay is determined by the distance along 
LOS, R x , giving by the first equation of Eq. Q, which is 
dependent of angles, i, rj and A, and distance, Rd- 

By Dermer (1999) the time corresponding to the radio 
peak is around 10 6 to 10 7 seconds, which means the distance 
of the hot spot to the core is, Rd = 10 16 to 10 lr cm. 

Consequently, the discrepancy in the time delay corre- 
sponding to the two hot spots, i.e., A2 and B2 in Fig. 2, can 
be given by putting two different precession phases, r?(ti) 
and r]{t2), corresponding to two different time, t\ and t^, 
respectively, into the first equation of Eq. (|4]), 

Ar _ R x (t2) - R x (ti) (n) 

c 

Therefore, by the parameters of 0119+041 and 0735+178 in 
the text, the time delay between two measured hot spots 
is less or around 1 day respectively, which is negligible in 
the fitting of Fig. 3 and Fig. 6. For GRS 1915+105 the time 



delay is ~ 3 hours which doesn't affect the fitting of Fig. 5 
either. 

If Rd exceeds lOOpc, then the time delay between two 
measured hot spots would be larger than 1 year. Under such 
circumstance, the hot spot, i.e., B2 in Fig. 2, would had been 
emerged at least 1 year before it was observed. In such case, 
increasing the jet precession speed, fi, by a few percent, the 
data can still be fitted. 

The two models can be tested simply by reanalyzing 
the data of 0119+041. If the previous features A and B 
move further away from the core when C is observed, then 
the ballistic model is supported, and the features should be 
produced by discrete jets. Otherwise, if any intermediate 
feature is found, i.e., between A and B, or B and C, then 
the non-ballistic model is favored, and the features should 
be powered by a continuous jet. 

The jet motion corresponding to a full period of pre- 
cession is expected on sources like 0119+041 and 0735+178. 
If this is observed then the non-ballistic model would be 
supported. 

The discrepancy of the inward and outward feature can 
also be tested. The apparent luminosity of an inward fea- 
ture increases gradually, and disappears in the core area 
(region 1), whereas that of an outward feature may decrease 
from a peak. And the core area appears moving towards the 
outward direction for certain angular distance, as shown in 
Fig. 7. 

The ballistic model is based on the assumption of a 
discrete jet, and can easily explain the linear jet motion. 
Comparatively, the non-ballistic model is based on the con- 
tinuous jet, which interprets both the linear and nonlinear 
jet motion. Therefore, it explains both the phenomena the 
ballistic model can explain, and the ones the ballistic model 
cannot. 

Testing the two models through their different predic- 
tions on both quasars and microquasars, and by both the 
data from the past and from the future, would enhance our 
understanding of the true mechanism of superluminal mo- 
tion, the interaction of jets with interstellar matter, and the 
common physics underlying quasars and microquasars. 
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